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Abstract 

Event-by-event measured arrival time distributions of Extensive Air Shower (EAS) 
muons are affected and distorted by various interrelated effects which originate from 
the time resolution of the timing detectors, from fluctuations of the reference time 
and the number (multiplicity) of detected muons spanning the arrival time distri- 
bution of the individual EAS events. The origin of these effects is discussed, and 
different correction procedures, which involve detailed simulations, are proposed 
and illustrated. The discussed distortions are relevant for relatively small observa- 
tion distances {R^ < 200 m) from the EAS core. Their significance decreases with 
increasing observation distance and increasing primary energies. Local arrival time 
distributions which refer to the observed arrival time of the first local muon prove 
to be less sensitive to the mass of the primary. This feature points to the necessity 
of arrival time measurements with additional information on the curvature of the 
EAS disk. 
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1 Introduction 



The pioneering studies of the EAS disk structure in the early 1950's by Bassi, 
Clark and Rossi [1] have been followed by increasing efforts (e.g. [2-8]) to de- 
tail the information, up to recent measurements with modern detector arrays 
[9-11] like KASCADE [12]. Especially for the EAS muon component [12-14] 
the basic interest arises from the question, how well does the temporal EAS 
structure reflects the longitudinal proflle and how well do muon arrival time 
distributions map the muon production heights, in particular when observed 
at sufficiently large distances from the shower axis, where the path-length ef- 
fects dominate [16]. Via these features muon arrival time distributions should 
carry some information on the mass of the EAS primaries, unless the intrinsic 
fluctuations of the muon generation processes and hmitations of the detector 
response do obscure the discrimination. Experimental studies with these as- 
pects are a current subject of the KASCADE experiment [17-20], taking use of 
the timing and muon detection facilities of the KASCADE Central Detector 
[21,22]. Experimental muon arrival time distributions have experienced vari- 
ous intriguing distortions which depend not only on the time resolution of the 
timing detectors, but also on the multiplicity, i.e. the number n of registered 
muons spanning the arrival time distribution of the single EAS event. Thus the 
shapes of the muon arrival time distribution are affected by the efficiency of 




Fig. 1. Characterisation of the EAS temporal structure by global and local arrival 
time parameters. 
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the detection system registering the fluctuating mimber of EAS muons. In this 
note we briefly discuss these multiphcity effects entangled with the response 
and the time resolution of the apparatus, and we suggest some procedures 
to reveal the basic time structure of the observed EAS. Fig. 1 indicates the 
used terminology (see also [17]). Arrival times of muons, registered by the 
timing detectors at a certain distance i?^ from the shower axis have to refer 
to a defined zero time. For the reference zero the arrival time Tc of the shower 
core in the detection plane could be used (global arrival times). The correction 
-^RijtanO / c (c - speed of hght) is apphed to muon arrival times in order to elim- 
inate the distorsions due to the shower inclination. However, there are often 
experimental difficulties to reconstruct this time-zero with sufficient accuracy. 
Hence alternatively local arrival times are considered which refer to the arrival 
time of the first muon, registered locally. In event-by-event observations the 
individual EAS relative arrival time distributions (global and local) can be 
characterised by mean values Ar^ean and by various quantiles like the median 
Ato.5, the first and the third quartiles Aro.25 and Aro.75, which describe dif- 
ferent features of the single distributions. The variation of the distributions 
of these quantiles, in particular of their mean values and the variances, with 
the distance from the shower centre, we call the EAS time profile. In the 
case of global time parameters these describe the curvature of the shower disk 
and the shower thickness, while local quantities characterise only the structure 
of the shower disk. The following discussion is mainly focused on the impli- 
cations of observations of local time quantities and is based on Monte Carlo 
simulations of the EAS development, using the code CORSIKA [23] with the 
hadronic interaction model QGSJET [24] as generator. 



2 Fluctuations of the arrival time of the first muon 

The question to which extent the first detected particle represents the arrival 
time of the EAS front has been discussed by Villiers et. al. [25]. 
Following these considerations it can be argued that the arrival time Ati of 
the first muon, relative to a fictitious zero time, representing the muon front, 
approximated by a sample with a large number n of muons, its expectation 
value and fluctuations depend on the particular value of the multiplicity n. 
Fig. 2 shows muon arrival time distributions accumulated from many showers 
with the expectation values (Ari(n)) of the arrival time of the first muon for 
subsamples of different multiplicities. 

With increasing n and in the case of infinite time resolution (cr = ns) of the 
timing detector (Ari(n)) approaches the fictitious arrival time of the shower 
front (which appears for small i?^ and infinite time resolution near Tc). In 
addition to the fiuctuations of the arrival times of the first muon due to the 
registered multiplicity (which in practical experimental cases involves also the 
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Fig. 2. Muon arrival time distributions ri(n) — Tc with the multiphcity dependence 
of (Ari(n)) and the influence of the time resolution. The shown distributions stem 
from the accumulation of 100 simulated Fe and 100 simulated proton induced EAS 
(zenith angle of incidence 0°) of 3-10^^ eV, observed in the range 70 m < < 80 m. 
In the present case of relatively small observation distances from the shower axis 
differences between Fe and proton induced showers would appear completely covered 
by the discussed fluctuations. 

influence of the detector response) there is the influence of the finite time 
resolution. It broadens the observed distributions and smears out the original 
asymmetry. As an example also the expectation for the (global) median (Ato.s) 
is shown. The averaged global median does not depend much on the time 
resolution. In contrast the local quantities underly the trends of (Ari(n)), i.e. 
they increase with the multiplicity and with the time resolution. 

These features are illustrated in Fig. 3 and lead to local shower profiles which 
are distinctly influenced by the detector qualities (time resolution and response 
for the muon multiplicity i.e. the detector area) as well as by the natural 
shower fluctuations. The influence is very pronounced in the central region of 
the shower with the largest multiplicities. As examples Fig. 4 displays local 
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Fig. 3. Multiplicity and time resolution dependence for the mean values of global 
and local time parameters of single muon arrival time distributions. 

time profiles (Atq.s) for proton induced showers observed with different time 
resolution. The profiles approach the ideal case at larger distances from the 
core (> approx. 250 m), where the time resolution loses obviously the influ- 
ence. It should be noted that local time profiles do exhibit nearly insignificant 
differences for different masses of the EAS primaries. This feature is contrast 
to observations of global time parameters [15]. 



3 How to account for the dependence on the registered muon mul- 
tipHcity 



In actual experimental observations of muon arrival time distributions like in 
the KASCADE experiment the profiles are derived from EAS events, includ- 
ing all different multiplicities (n > n^,), as registered by the timing detectors. 
Due to the lateral distribution of the EAS muon component the average (reg- 
istered) multiplicities in an EAS are depending on R^. They depend also on 
the type and energy of the primary, and even if the energy is approximately 
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Fig. 4. Local time profiles of proton induced EAS for different primary energies and 
different time resolution, virtually registered with a detector eye of 400 detectors 
of 0.5 (approximately like the timing facility of the KASCADE central detector 
[12]). 

specified, the observed time profiles originate from a superposition of various 
multiplicities, varying with the distance from the shower core. That feature 
leads to interference effects, distorting the predicted quasi-parabolic shape 
of the time profile. Fig. 5, which displays measured EAS time profiles [26], 
exhibits this effect. 

In principle, it would be desirable to extract the shower profile for each multi- 
plicity separately. However, this attempt would meet serious problems of the 
statistical accuracy of the results. In addition such a procedure would lead to 
a less transparent selection of shower events. The main problem, however is 
that the appearance of the distortions depends on the response qualities of the 
particular detector arrangement, so that measurements by different detector 
arrays are hardly directly comparable. There are various ways to approach 
a representative result about the EAS time structure from measurements of 
local quantities. 

• The standard concept to compare experimental results with theoretical pre- 
dictions are procedures simulating the experimental conditions and response 
folded with the predictions and comparing the resulting distributions with 

the measured results. This will give an impression about agreement or dis- 
agreement and on the undistorted time profiles, but does not immediately 
enable the comparison with other experiments of different (often unknown) 
quality and conditions. 
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Fig. 5. The observed variation of the mean values of the median, first and third quar- 
tile distributions with the distance from the shower axis (from [26] ) , extracted with- 
out multipUcity correction for an EAS sample within a particular range of log^g 
(indicating [27] the primary energy range of approx. 1.6 • 10^^ eV to 3.2 • 10^^ eV). 
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The measured time parameters, deduced for each event from the single muon 
arrival time distributions registered with varying multiplicity, get scaled to 
a chosen reference value of the multiplicity by a correction according to pre- 
dictions by simulation calculations (indicated for the case of the mean value 
Axmean by a Calibration curve in Fig. 3). The appearance of time profiles de- 
pends on the choice of the reference multiplicity. The correction procedure 
needs detailed simulation calculations of muon arrival time distributions of 
the kind shown in Fig. 3. The procedure has been successfully applied in 
recent KASCADE experiments [17,19]. 

The observation that for global time quantities the influences of the multi- 
plicity and the time resolution are less pronounced, suggests to relate the 
muon arrival times to the arrival time of the shower centre by simulating the 
time difference between the arrival time ti of the local first muon and the 
arrival Tc of the EAS core. In this way [20] the local quantities are trans- 
formed into pseudo-global time parameters, which display the EAS time 
structure rather realistically, but invoke EAS simulations, specified in de- 
tail (see also ref. [11]). While only the shape of the arrival time distributions 
enters in the above multiplicity calibration procedure, the transformation 
to pseudo-global quantities stresses also the absolute time difference ri — Tc 
between the arrival of the first muon and of the shower core in the plane 
perpendicular to the axis. 
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Fig. 6. Local EAS shower profiles ((Aro.5)), corrected for the multiplicity depen- 
dence, as compared with predictions of EAS simulations using the Monte Carlo code 
CORSIKA [23] and taking into account the time resolution, for different EAS muon 
sizes (corresponding to an energy range from 6 • 10^^ eV to 6.3 • lO-*^^ eV). 

Finally Fig. 6 presents a result of an experimental investigation [19] of EAS 
time profiles using the KASCADE detector. For different ranges of the trun- 
cated muon number Nj[ (used as approximate energy estimator [27]) and 
with a consistent correction for the multiplicity dependence, the experimen- 
tal results are compared with CORSIKA [23] simulations, adopting a mass 
composition p:0:Fe = 1:1:1. 



4 Concluding remcirks 



The shape of measured muon arrival time distributions, which refer to the 
arrival time Ati of the first locally registered muon, experiences some dis- 
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tortions which arise from the locally varying superposition of different muon 
multiplicities (spanning the observed individual event distribution). The effect 
of multiplicity fluctuations is entangled and amplified in an intricate way with 
the response and time resolution of the detector setup. Due to the larger inten- 
sity of the electromagnetic EAS component and the smaller fluctuations of the 
local arrival time of the first charged particle, the effects are less pronounced 
in practical cases of studies of the time structure of the charged particle com- 
ponent [12]. In principle, however, the dependence of (Ari) from the observed 
multiplicity is also present in simulation studies with an ideal time resolution, 
though mostly covered by the assumption of large collection areas of the ideal 
detectors. There are various procedures, which remedy the distortions and 
reveal the basic experimental results of the structure of the EAS disk for a 
reasonable comparison with theoretical predictions or other experiments. The 
described procedures are weakly dependent on the used interaction model, 
governing the EAS development, and on the adopted mass composition. They 
invoke explicitely EAS Monte Carlo simulations. The described distortions of 
the observed muon arrival time distributions lose their significance with in- 
creasing distance from the shower axis ( > 250 m) and with increasing 
energy. Nevertheless local time distributions (referring to the arrival time of 
the first local muon) turn out to be rather insensitive to the mass of primary 
cosmic rays. This is in contrast of global arrival times [12-15] which refer e.g. 
to the arrival time of the shower core and measure additionally the curvature 
of the shower disk. 
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